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Accurate expressions, valid in experimentally relevant regimes, are presented for the effect of long- 
ranged Coulomb interaction on the low energy properties (momentum distribution function, density 
of states, electron spectral function, and 4k F correlation function) of one-dimensional electron sys- 
tems. The importance of plasmon dispersion (as opposed to exponent) effects in the spectral function 
is demonstrated. 
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The low energy behavior of one-dimensional (ID) elec- 
tron systems is known not to be consistent with Fermi 
liquid theory M. However, the theoretically well estab- 
lished H and extensively studied J| Luttinger liquid 
model of one-dimensional physics pj is, strictly speak- 
ing, not applicable to electronically conducting one- 
dimensional systems such as quantum wires (QWR) ||, 
carbon nanotubes organic conductors and doped 
chain or ladder compounds because the electrons in 
these compounds interact via the Coulomb force, which 
is long-ranged, whereas the standard Luttinger model as- 
sumes short-ranged interaction. The long range of the 
Coulomb interaction leads to a scale dependence of the 
Luttinger exponents and velocities which have been 
studied by several authors |&H on the assumption that 
it is well approximated by its leading (In 1 / 2 ) divergence. 
As we show, this approximation is not accurate in any 
physically relevant regime. One exception is a very inter- 
esting recent renormalization group treatment jll ] which 
found an effective exponent very similar to ours but did 
not discuss the implications for physical quantities. Some 
numerical results have also appeared [ fl2"| , but a general 
understanding of the experimental implications of the 
Coulomb interaction is lacking. 

In the present paper we use direct analytical and nu- 
merical evaluation of the relevant bosonization expres- 
sions to determine the momentum distribution function, 
tunneling density of states, and spectral function for 
ID electron systems interacting via the physically rel- 
evant Coulomb interaction at zero temperature. We de- 
fine an important but previously overlooked energy scale, 
present an accurate expression for the scale dependent 
exponent, show how the scale dependent velocity affects 
the spectral function, and qualitatively discuss the 4k F 
correlation function. Our results should apply directly to 
ID QWRs § and nanotubes §. 

We consider a ID electron system with a noninteract- 
ing dispersion s p which we linearize near the Fermi point, 
defining a bare velocity v F . We here assume the only im- 
portant interaction is the Coulomb interaction in the for- 
ward scattering channel, and neglect Umklapp scattering 
and other interactions. This is a good approximation for 



QWR and nanotube systems. (For the organic and doped 
spin chain materials a modification, discussed below Eq. 
(2), is needed.) The Hamiltonian is (here we do not write 
the spin index explicitly), 



H = ^2v F (p - rp F )cl. 

r,p 

+7 Vc ^ (pr(q)pr(-q) + pr(q)p-r(-q)) , (i) 
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where c\ p is the electron creation operator and p r (q) 
is the density operator describing density fluctuations 
at momentum q and branch r = ±1 for the right(lcft) 
movers. For ID systems V c (q) — ► ln(l/g) as q — > 0, 
and becomes 1/q for q larger than some scale qo se t by 
the geometry and the wave function size. A reasonable 
approximate form, which we will use in our subsequent 
analysis, is 
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where Vq is a dimensionless measure of the interaction 
strength and q^ 1 is the length scale parameter. Vq and 
qo are system dependent factors. For a cylindrical quan- 
tum wire of radius a, Vq — 4e 2 /ireovp and qo ~ 2.5/a, 
where e is the electron charge, and Eq is the background 
dielectric constant, about 10 for GaAs. These values give 
the correct long wavelength limit and are within 10% of 
the correct 1/q coefficient at large momentum. In car- 
bon nanotubes Vq is of the same form as in QWR but 
Eq ~ 1.4 |J and qo ~ 2.97/i?, where R is the radius 
of the tube. For organics or doped spin chains, addi- 
tional short-ranged exchange interactions may be impor- 
tant. The usual arguments [Q show that these interac- 
tions lead, at low energies, to an additive constant term 
in V c (q). 

Eq. (1) may be bosonized as usual j2|||; the charge 
excitations are plasmons with dispersion uj q — qv q , and 
velocity v q = v F \/i + 2V c (q)/irv F is 
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(Note that we have lim^ooU^ ~ qvp + Voqo/2 7^ qvp 
for Coulomb interaction). The electron Green function 
G r (x,t) = (il) r (xt)ipl(00)) is 

G r ( X ,t) = ^n eirk ' XieXp[ -^ X ' t)] , (4) 
v ; £^0 2tt x-rvpt + ie v ' 

where the phase function & r (x,t) is 

1 f°° dp 
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The exponent parameter 9 p is defined by 



rvpt) gip(x — rv p t) 



2 sinh 2 (1 - cos(px)e' lrpv " t ) . (5) 
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where q s = q^e 1 ^ and the last approximation is good 
at long wavelengths, q <C qo- 

We now use Eqs. (4)-(6) to study electronic quantities. 
We begin with the momentum distribution function 



n r {8p) = — 



dx- 



-idpx 



■exp[-$ r (x,0)], (7) 



where Sp = p — rkp. In a noninteracting Fermi gas, 
n r (Sp) — 0(—rdp). For a short-ranged LL, the generally 
accepted result |l3| is that in the vicinity of the Fermi 
momentum, 0.5 — n r (Sp) ~ sgn(rSp) x {Ci\5p\ + C2\Sp\' 1 } 
with 7 a LL exponent and C\ and Ci two constants. 
The first term is the non-critical background coming from 
high energies, while the second (critical) term comes from 
low energies where LL physics is important. For the long- 
ranged interacting model we now consider, attention to 
the singularity structure of the noninteracting electron 
Green function leads to (let 8p > and r = +1) 

1 1 f°° dx 

n(Sp) = / — sin (Spx) exp[— $(x, 0)1 

2 n J x 

1 {5v\ lq(Sp) 

= -+ C[5p + C'J — ) + higher orders, (8) 
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where again the non-singular C[ term is from the inte- 
gration over small x, while the singular C' 2 term comes 
from integration over large x and is a weak function 
of ln 1 ^ 2 (l/<5p). The scale dependent exponent r ) q {Sp) is 
found to be 
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Fig. 1 shows results obtained by numerically eval- 
uating Eq (8) for typical QWR parameters. An en- 
hanced curvature near the Fermi momentum is evident. 
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FIG. 1. Calculated momentum distribution function, 
n(p), with respect to momentum p — kp for a realistic 
QWR system of a — 70 nm. Solid, dashed and dot- 
ted lines are results for three different interaction strengths 
Vo = 1.21,2.42,4.84 respectively, where Vo = 1.21 is for elec- 
tron density 0.65 x 10 6 cm' 1 and eo = 12.7 [5]. Inset: the 
effective exponent (around fcj?) obtained by taking the loga- 
rithmic derivative of the numerical n(p) for \p — &f < 0.2fcF- 



The inset to Fig. 1 shows the logarithmic derivative 
a q (p) = dlog\n(p) — l/2|/dlog(p), which shows that for 
small Sp the behavior may be described in terms of a 
slowly changing effective exponent. We note a q (p) is al- 
ways less than 1, because when the scale dependent ex- 
ponent 7 g (p) of Eq. (9) is greater than 1, the background 
term dominates. 

We now turn to the tunneling density of states, 



N(lo) 



1 f°° 

—J2 dt e ' ui [ G -(°> *) + G -(°> -*)] . ( 10 ) 



for to measured from the chemical potential. We first 
show that N(uj) vanishes faster than any power of uj as 
uj — > 0. We observe that if V c (p) 7^ 0, G vanishes faster 
than any power of t as t — ► 00 M . Therefore the integral 
obtained by taking any number of ^-derivatives of N(u>) 
is absolutely convergent at long times, and may be evalu- 
ated straightforwardly by contour methods even at u = 
|l4j| . We further note that G>(0, t) has no singularities 
in the lower (upper) half plane for r = +1(— 1); thus 
by deforming the contours appropriately we find that 
d n N(u;)/dui n \ u , = o = for any n. This argument does 
not apply to n(p) because of the different analytic struc- 
ture of the x-dependence. Thus the non-critical contri- 
butions which obscured the behavior of n(p) do not occur 
in N(lu). By evaluating Eq. (10) we obtain 



N(lo) oc ( — 
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where the scale dependent density of states exponent 
7^M is 
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FIG. 2. Calculated density of states, N(ui), with respect 
to energy u> for the same system as Fig. 1. Different line 
styles represent different interaction strengths as indicated. 
Inset is the effective exponent, a u obtained by taking the 
logarithmic derivative of N(ui). The numerically calculated 
curves are well fitted by the analytical expression (dash-dot 
fines) of the exponent from Eq. (13) at the corresponding 
Vb's and ui s = 20g s «FV / Vo- The stars are the first order term 
of Eq. (13) only, for Vb = 1.21, showing that the widely used 
leading logarithm approximation leads to factor of two errors. 
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the same form as that of Eq. (9) with q s replaced by a 
characteristic energy scale ui s . From Eqs. (3) and (6) 
we expect ui s = Aq s VF \/Vo with the numerical constant 
A determined by subleading corrections to the asymp- 
totic analysis of Eq. (5). A may in principle have a 
weak scale and system-parameter dependence, but our 
numerical results show that for a wide range of energies 
(1CT 3 < lu/E f < 0.1) and interactions (1 < Vq < 5) it is 
very well approximated by the constant value A = 20. 
Fig. 2 shows the results of a numerical calculation 
of N(u>) from Eq. (10) for three different interaction 
strengths; the inset compares the numerically calculated 
effective exponent, a u (u) = d\og(N(oj))/d\og(uj), with 
the analytical result obtained from Eqs. (11) and (12): 
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One sees that the fit is very good (the small differences 
appearing at ui/Ef ~ 0.01 arise from noise in the numer- 
ical calculation). 

The two crucial energy scales defined by N{uj) are 
uj 8 , and at which a u (uj*) — 1, corresponding to 
~ uj s e~ 34 / v ° <C ui s . In the high energy region, 
uj > oj s , one has essentially non-interacting behavior. 
For to* < lo < u> 8 , one has a LL with a scale depen- 
dent exponent. For uj < uj*, > 1 and N(u>) is concave 
upwards at small uj, suggesting a "pseudo-gap" in the 



electronic density of states. For most real QWR systems, 
Vb is about 1-5 depending on Eq and v F , and thus uj* is 
typically many orders of magnitude smaller than w,. In 
our calculation, using QWR parameters from Rcf. If, we 
have uj s ~ 100 meV and a->* ~ 10~ 4 meV. For extremely 
small uj < w», Eq. (13) gives a w ~ ^ ln 1/2 (uj s /uj), 
an approximate form used earlier in the literature ||,[l0| • 
However, as seen from the inset of Fig. 2, the leading 
logarithmic divergence is so weak that in all physically 
relevant regimes the other two terms in Eq. (12) are 
needed for quantitative accuracy. On the other hand, 
the constant (scale independent) exponent used in Ref. 
H for nanotubes is also not an adequate approximation 
for small energy region (uj < 0.05Ep) either. We there- 
fore propose Eqs. (9) and (12) as widely applicable fit- 
ting formulae for the effective exponents in the Coulomb 
Luttinger liquid. 

The scale dependent exponent also appears in the sin- 
gle particle spectral function, p(q,uj) = (l/2ii)[G{q,uj) + 
G(—q, —oj)], however the scale dependent velocity in Eq. 
(5) is more important. To introduce our results, we 
briefly summarize known results for a short-ranged re- 
pulsive interaction in the spinless LL model |l3|] . At fixed 
q, one defines three w-ranges: (i) p(q,uj) = for \u>\ < 
uj q (energy-momentum conservation), and (ii) power-law 
singularities as |o>| — * and (hi) an exponential decay 
at scales larger than the Luttinger cut-off. For the long- 
ranged Coulomb interaction, p(q, uj) = for \uj\ < ui q due 
to the energy-momentum conservation, but the behavior 
in both regions (ii) and (hi) are strongly modified. For 
|u;| > uj s (region (hi)), p(q,oj) ~ exp[— |a>|/.E c (w)] with a 
scale dependent cut-off 



„ , \ qoVpVo , f a\w\ 
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because of the slow (l/q) decay of the Coulomb interac- 
tion in the large momentum region (Eq. (5)). 

Near threshold (uj q < \w\ <C uj s ) there are two effects: 
the scale dependence of the effective Luttinger exponent 
and the curvature of the plasmon dispersion, which pre- 
vents the different boson modes from adding coherently. 
Thus as one decreases uj towards uj q (consider uj > 
part only) one obtains first a divergence 

5u p„(Su)-i ( here 

5u> = uj — uj q ). This divergence is cut off by curva- 
ture effects at a scale uj c (q) = Max p<g (w p —puj q /q) w 

[1/ 4)qvF\/Vo ln _1 ^ 2 (<Zs/'?), the difference between the ex- 
act dispersion and a linear approximation. We find that 
for q larger than g» ~ Aq s e~ 75 ' v ° (at which 7 a ;(ij- l c (g*)) = 
1) the curvature effect is more important in cutting off 
the divergence, whereas for q < q* the effective exponent 
is more important. As Suj — > + the spectral function 
decreases rapidly, ultimately vanishing faster than any 
power of Suj due to the increase of the effective exponent. 
Thus the generic behavior is a spectral function which in- 
creases rapidly as uj is increased above threshold u> q , goes 
through a maximum at uj pea k = uj q + Auj with Aw set by 
the larger of w* and u) c (q), and then decreases exponen- 
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FIG. 3. Calculated electron spectral function, p{q,u)), for 
different momenta, q, as indicated in the figure. Solid lines 
are from the Coulomb interacting system (parameters are the 
same as Fig. 1 with Vo — 1.21) while dotted lines are from 
short-ranged interacting system with the approximate effec- 
tive exponent a = 0.2. The two triangles in (b) and (c) 
indicate the threshold, u = ±uo q . The ripple of the spectral 
function curves is numerical error. 

tially with a scale dependent cut-off E c (uj) for w > u> pea k- 
The suppressed spectral weight in the near threshold re- 
gion is compensated by the slower decay at high energies, 
preserving the sum rule J p(q,tu)duj = 1. In Fig. 3 we 
show the results of direct numerical evaluation of the elec- 
tron spectral function for the Coulomb Luttinger liquid 
(solid lines) and for a short-ranged-interaction (regular) 
Luttinger liquid with exponent a = 0.2, approximately 
equal to the effective exponent of the Coulomb case at 
u> = Q.SEp (dashed lines). Note that the partition theory 
techniques used in [ fLSj to simplify the evaluation for the 
short-ranged case do not work in the Coulomb case. The 
shift of the peak away from the threshold is evident. 

Finally, we briefly discuss the "Wigner crystal" cor- 
relation. Schultz H observed that at long enough 
length scales the logarithm arising from the long-ranged 
Coulomb interaction causes the 4fc^ component of the 
density-density correlation to decay more slowly than 
x -4if p an( j a i g0 more sl ow ly than the 2kp component, 
leaving a state best interpreted as Wigner crystal. Using 
the notations of this paper, we obtain for the term 
in the structure factor, 



/ n l-8Vo- 1/2 ln- 1/2 (W5p) 

S <**( S P) ~ , (15) 

where dp = \p— 4fcjr|. Therefore we expect to see the Akp 
divergence when ^/Vo In 1 / 2 (q s /5p) > 8 or Sp < q s e~ 64: / Vo , 
or in term of temperature at T < T w , x = uj s e~ 64: / Vo , 
which is sensitive to the electron density and experimen- 
tal geometry, but is in general far too small to be exper- 
imentally relevant, and is also much less than the scale 
at which N(u>) develops a pseudo-gap. 

In conclusion, we have presented a systematic theo- 
retical analysis of the low energy properties of electron 
systems subject to long-ranged Coulomb interactions, in- 
cluding a reliable estimate of the scale dependent Lut- 
tinger parameter and apparently the first calculation of 
Coulomb effects on the spectral function, and values for 
the (unfortunately extremely low) scales at which the 
divergent behavior associated with the Coulomb interac- 
tion becomes manifest. 
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